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We analyze theoretically a possibility of the superradiant emission (SR) in GaN-based quantum 
well (QW) laser heterostructures. Our model is based on travelling wave Maxwell-Bloch equations 
predicts building up of macroscopic coherences in the ensemble of carriers in the QWs. We show that 
SR is covered by the Ginzburg-Landau equation for phase transition to coherent matter state. The 
generation of superradiant pulses as short as 500 fs at peak powers of over 200 W has been predicted 
for InGaN/GaN heterostructure QWs with the peak emission in the blue/ violet wavelength range. 
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A spontaneous build-up of macroscopic coherences in 
solids, or in other words, a phase transition to a single 
macroscopic quantum state of matter, has always been 
a fascinating subject for scientific and engineering com- 
munities. Recent examples of quantum coherent phe- 
nomena include the Bose-Einstein condensation (BEC) 
of atoms, which stipulated works on laser-cooled atomic 
clocks with unprecedented fractional frequency stability, 
and wide spreading use of superconducting magnets in 
which a coil undergoes the BCS (Bardeen, Cooper, and 
SchriefTer) transition. The search for spontaneous build- 
up of macroscopic coherences in semiconductor microcav- 
ities has resulted in intensive research towards a polariton 
laser. Interestingly, group III - nitride semiconductors 
are capable of reaching macroscopic coherences at room 
temperature conditions. Thus due to large reduced exci- 
ton mass, the InGaN/GaN quantum well (QW) exciton 
binding energy is higher as compared to GaAs/AlGaAs 
QWs or other conventional III-V counterparts. Further- 
more, due to a stronger exciton-photon Rabi coupling, 
InGaN/GaN microcavities offer possibility of polariton 
lasing at room temperature conditions [1 j. High critical 
temperature of BEC transition is conditioned by low in- 
plane effective mass of polaritons, so as the critical den- 
sity is reached before destruction of excitons. Such 2D 
BEC is a transient dynamic state, which can nevertheless 
be analyzed using a stationary Schrodinger equation [2]. 

The superradiance (SR) in semiconductor edge emit- 
ters can be considered as another example of spontaneous 
macroscopic coherences in solids [3] . The cooperative ra- 
diative recombination in an ensemble of quantum oscilla- 
tors (e.g. atoms or molecules) has been predicted before 
the invention of lasers [4 j. Since that, it has been ex- 
tensively studied both theoretically and experimentally 
[5j [6] . The characteristic features of the SR emission are 
the temporal and spatial coherence, highly anisotropic 
emission pattern, quadratic dependence of pulse intensity 
on the number of excited atoms I oc n 2 , afterpulse ringing 



attributed to Rabi-type oscillations. The SR pulse du- 
ration decreases with the number of emitters r sr oc 1/n, 
while the pulse energy is proportional to the ensemble 
population (Ir sr oc n), in agreement with the energy con- 
servation considerations. On the other side, the sponta- 
neous nature of the transition to the transient macroscop- 
ically coherent state is responsible for large fluctuations 
in the shape, duration and amplitude of SR pulses due 
to quantum- mechanical uncertainties [7J |8] . 

In semiconductors, a hypothesis has been drawn that 
the SR is assisted by formation of a transient coupled 
electron-hole (e-h) pair state mediated by photons, much 
like as phonon mediation in a superconducting BCS state 
[9j [10] . According to that the e-h system undergoes the 
second order non-equilibrium phase transition when the 
coherent e-h BSC-like state is building up during SR 
pulse emission. From practical point of view, SR has 
been considered as one of the promising approaches to 
generate high-power femtosecond optical pulses exceed- 
ing by a few orders of magnitude the power in lasing or 
amplified spontaneous emission regimes. 

In this letter, inspired by remarkable features of 
the group III - nitride microcavities in reaching room 
temperature BEC macroscopic coherences, we stud- 
ied numerically the Dicke superradiance in an edge- 
emitting ridge- waveguide multiple-section cavity with In- 
GaN/GaN QWs (see Fig [I]). We developed a numerical 
model based on a semi-classical description of macro- 
scopic coherences (polarization) in the ensemble of elec- 
trons and holes in the QWs. In distinguishing from pre- 
vious theoretical treatment developed for the bulk and 
QW GaAs-based materials [11 , our new numerical model 
takes correctly the carrier populations and coherences 
into account in function of the pump current and ab- 
sorber bias. To support the hypothesis of BCS-like con- 
densation, we show that the travelling wave Maxwell- 
Bloch equations governing the SR emission can be con- 
verted to the Ginzburg-Landau equation for the second- 
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order BCS phase transition. This enables us to define the 
critical density in the system, the order parameter and 
the coherence time (and length). 

As a model system we utilize a tandem cavity config- 
uration consisting of two separately contacted sections 
(Fig. [I]). The shorter negatively biased section is used as 
a saturable absorber. The longer section is pumped with 
short current pulses and provides the optical gain. 

The evolution of the amplitudes of the forward 
and backward (A-) travelling waves is defined by the 
cavity loss and macroscopic medium polarization 



± v Q 



-v g aiA±, 



(1) 



dt y dz 2 V T 2 
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The normalization used for coherences (the ensemble av- 
erage for the off-diagonal density matrix elements) in- 
duced between the electrons and holes and associated 
with the forward and backward waves reads: 
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In convention (J3|, the macroscopic variables P± measure 
the order parameter of the system. The usual bound- 
ary conditions are applied at left (z=0) and right (z=L) 
cavity mirrors 

A+(0, t)=VRA_ (0, t), A_ (L, t)=VRA+(L, t). (4) 

Separate equations for the gain and absorber sections 
take into account that (i) the carrier lifetime in absorber 
is much shorter as compared to the gain; (ii) pumping 
by carrier injection is applied to the gain section only, 
(hi) Quantum Confined Stark Effect (QCSE) absorption 
tuning occurs in the absorber section. The dynamics of 
carrier populations and coherences is described as fol- 
lows: 
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where the carrier diffusion D is taken into account; a is 
the differential absorption to gain ratio. The Langevin 
force term due to polarization (thermal) noise A± triggers 
the spontaneous build-up of superradiance pulse. 

Figure [2] displays the numerical solution of the above 
equations for the SR emission in the cavity of the overall 
length L=600/im and relative absorber length 20%. The 
absorber section is located at the left facet of the struc- 
ture (at z=0) and the gain section is situated at the right 
facet. 
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FIG. 1: (a) Schematic of two-section laser structure for SR 
emission generation and (b) output SR pulses at the gain sec- 
tion facet (red curve) and absorber section facet (blue curve) . 



Initially, a high reversed bias is applied to the absorber 
leading to depletion of carriers in the QWs (blue color 
in the scale of Fig{2|a)). It prevents the structure from 
lasing and enables accumulation of carriers in the gain 
section. The gain section QWs are pumped to a high 
level, so as in this particular case, the carrier population 
is about 15 times above the transparency n t (red color 
in Figj2|a). The entire process starts at t=0 from the 
spontaneous polarization noise. The rate of spontaneous 
polarization noise into the cavity mode ~^rV#o^2A± 
corresponds to 27 photons per cavity round trip. This 
cannot be seen in the scale of Fig(2| The noise sources 
are uniformly distributed along the cavity but the macro- 
scopic coherences are building up at the edges of the 
gain section, after about half of the cavity roundtrip 
time (T cav =14ps). Because of the asymmetric cavity, the 
shape and emission time of the output SR pulses from 
the left and right cavity facets are different (FigJTJb)) 
and the emitted SR pulses are not related by reflections 
at the cavity facets. 

The build-up of the macroscopic polarization is fol- 
lowed by Rabi oscillations at a frequency dependent on 
the pumping rate . They are clearly seen in the wave am- 
plitudes (Figsgd) and (e)). The carrier density drops 
abruptly at the same time when the polarization rises 
indicating that almost all carriers contribute to the field 
(Figs|2ja) and (c)). When the travelling backward pulse 
hits the absorber section at z=0.2L, the pulse intensity 
is sufficiently high to saturate absorber (the blue step in 
the region 0<z<120/im at t=7ps in Figj2^a)). From this 
time, there is no absorption and gain in the cavity. The 
backward travelling SR pulse reflected of the left cavity 
facet can now freely travel through the structure yield- 
ing emission of the secondary pulse at the right cavity 
facet (Figjljb), red curve). In this particular realization 
of SR emission, the pulse emitted form the absorber sec- 
tion side is shorter due to sharpening of its rising front 
in saturating absorber. 

The build-up of macroscopic polarization has a 
stochastic nature, which results in differences between in- 
dividual realizations of SR pulses. The peak power and 
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FIG. 2: Spatiotemporal dynamics of carrier populations (a), 
macroscopic polarization (carrier coherences) in InGaN/GaN 
QWs associated with the forward (b) and backward (c) waves, 
field amplitudes of the forward (d) and backward (e) waves. 
The axes T and Z are normalized on cavity roundtrip time 
and cavity length, respectively. The polarization and the field 
amplitudes are normalized on saturation parameter. The cur- 
rent amplitude I = 211 mA. 



pulsewidth of SR pulses vary within the range of 50-200 
W and 0.5-5 ps, respectively, depending on the driving 
conditions (Figj3|. Varying the overall cavity length and 
the relative length of absorber, we find that the inten- 
sity ratio and the width of the forward and backward SR 
pulses can be effectively altered. The common feature of 
all realizations is that the polarization and carrier excess 
above the transparency vanish after a few oscillations so 
as the the optical pulses travel freely without gain or ab- 
sorption. 

In Figj2](e), the peak intensity of the pulse corresponds 
to photon density in the cavity of 9-10 18 cm -3 . The den- 
sity of the associated coherence excitations P± at the 
polarization peak maxima is of 2.8T0 21 cm -3 , which is 
much higher than the initial carrier density no=2.4-10 20 
cm -3 . A hypothesis has been previously stated [9l [10] 
that the SR in a semiconductor is mediated by a forma- 
tion of the transient BCS-like state of e-h pairs. 

In order to interpret the SR emission in terms of BCS- 
like transition one should transform Eqs.Q and ([5| to 
the form of the Ginzburg-Landau equation (GLE) for 
the second-order BCS phase transition. The difficulty 
is caused by the fact the GLE is essentially a steady- 
state equation for an isolated system, while the SR is a 
transient process. To overcome this difficulty we change 
variables to the internal coordinates of the pulses. For 
simplicity we consider only the forward traveling pulse, 
introducing new coordinate system (=t — z/v g and z. 
The phase transition under the question is related to the 
change of the effective phase relaxation time T| ff , which 
is comparable with the lifetime of photons (l/i^c^=3ps). 
At the same time, the carrier relaxation can be neglected 
at the time scale of the pulse width: 
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The analytical solution is obtained for a model with local- 
ized parameters. The substitution 1/L* — » d/dz yields 
A+=P+r^go/T 2 /v g (ai+2/L*) with L* being some ef- 
fective length (see the discussion below), which allows us 
to exclude the wave amplitude: 
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The last equation recovers the Bloch vector conservation 
(n(C)— n cr ) 2 +P + (C) 2 = (^o— n cr ) 2 with no = n{— oo) be- 
ing the initial carrier density in the system. Finally, tak- 
ing the second derivative in the first equation ^ one 
obtains the master equation: 

d 2 P + T 2 glL* 2 



d( 2 T 2 v 2 g (2+aiL*) 2 



(n -n cr y-2P^ P + =0. (8) 
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This equation has a canonical form of the Ginzburg- 
Landau equation (GLE), in which P + is the order param- 
eter of the system and the internal coordinate (=t—z/v g 
plays the role of spatial variable. For the carrier densi- 
ties above the critical n cr (see Eq.([7])), its steady state 
solution defines the condensate fraction (BCS-like). The 
analytical similarities between Eq.([8| and GLE, allows 
one to assume that for n>n cr , the evolution of the the 
order parameter is defined by hyperbolic secant function. 
Substituting P + =Po secn (C/ r p)? we obtain the coherence 
time r c =21n(l+V / 2)r p (the FWHM of the SR pulse), and 
the peak photon number in the cavity (7V + =A^_(0)) 

T P=-r^h ^> N += o ^ ( n o-n cr ) (9) 

Tgo{n -n cr ) gor* 

The obtained pulse width 2\n(l-\-^/2)r p and peak power 
oc Tiujv g Nj r exhibit the initially expected behavior for the 
SR pulse known from the studies of SR in a gas medium. 
We thus shown that (i) the SR in semiconductor is the 
BCS-like phase transition and follows the GLE ^ and 
(ii) the order parameter of the system showing the abrupt 
growth is the medium polarization P±, which is defined 
by coherences (off-diagonal terms) in the density matrix. 
This indicates that the wavefucntion of the condensate 
quasi particles is the composed of products of the e and 
h wavefucntions. The proposed theory of BCS-like con- 
densation during the SR emission is different from all 
previous treatments. Although it does not involve the 
system temperature explicitly, the thermodynamic tem- 
perature can be introduce to characterize the the degree 
of population inversion, that is the initial quasi equilib- 
rium state of carriers of density no- 

The SR pulse features ([9| are in perfect agrement with 
the predictions of the Dicke theory if the coherence length 
L c =T c v g is longer than the cavity length L. In that case, 
the characteristic length of the condensate fraction L* 
is defined by the sample size L*=L<L C . However with 
increasing carrier density no, the coherence length L c be- 
comes shorter than the cavity length L and the size of 
the condensate fraction is set by the coherence length 
L*=L C <L. As can be seen from Fig. [2](b) and (c) show- 
ing the order parameter P±, this is the typical situation 
in the semiconductor laser cavity. Thus for n > n cr 2 [32] , 
Eq.([9| leads to a quadratic equation with respect to r p . 
Its asymptotic behavior for n ^> n cr2 reads 



N+ cc (n -n £ ) 3/2 
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p y Tg (n-n t y 
It is interesting to compare the predictions of our sim- 
plified analytical model with the results of our numer- 
ical experiment in Fig(2j The critical carrier density of 
n cr =2.52n t is achieved at 4.2 kA/cm 2 pump current den- 
sity (n=r n J/e q d) so as the condensation condition n>n cr 
is fulfilled. However at n cr2 = 2.57n t (at J=4.3kA/cm 2 ) 
the coherence length reduces down to the size of the cav- 
ity. Thus only in the narrow range n cr < n < n cr2 , the 
SR pulse peak power exhibits Dicke's quadratic growth 
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FIG. 3: SR pulse peak power (a), pulse width (b). 



([9|. At higher injection levels, it is switching to the de- 



pendence with the asymptotic (10). This feature is per- 



fectly reproduced by the numerical model in Fig. [3| In 
conditions of Figj2j the estimated pulse width 1.76r p is 
600fs. 

The L-I and pulse width curves in Figj3] exhibit large 
amplitude and timing instabilities of the output SR 
pulses. The SR is particularly interesting quantum op- 
tics phenomenon because quantum fluctuations manifest 
themselves macroscopically in the time and energy do- 
mains. The predicted macroscopically large fluctuations 
are one of the characteristic features of a condensed state 

In summary, we presented the numerical model of 
Dicke SR in semiconductor laser cavities and shown that 
it corresponds to the analytical Ginzburg-Landau equa- 
tion for BCS phase transition. We have demonstrated 
that the theoretically predicted performance of SR pulses 
in InGaN/GaN QW laser heterostructures exceeds sub- 
stantially that of the mode-locked or Q-switched GaN- 
based lasers [13] . 
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